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Abstract 

We complete the effective potential calculation of the two-loop, top/bottom Yukawa cor- 
rections to the Higgs boson masses in the Minimal Supersymmetric Standard Model, by 
computing the 0(a^ + ottocb + o^) contributions for arbitrary values of the bottom Yukawa 
coupling. We also compute the corrections to the minimization conditions of the effective 
potential at the same perturbative order. Our results extend the existing O(af) calcula- 
tion, and are relevant in regions of the parameter space corresponding to tan/3 3> 1. We 
extend to the Yukawa corrections a convenient renormalization scheme, previously proposed 
for the C(af,a s ) corrections, that avoids unphysically large threshold effects associated with 
the bottom mass and absorbs the bulk of the corrections into the one-loop expression. For 
large values of tan j3, the new contributions can account for a variation of several GeV in the 
lightest Higgs boson mass. 
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1 Introduction 



One of main features of the Minimal Supersymmetric Standard Model (MSSM) [T] is the pre- 
diction of the existence of at least one light Higgs boson 2j. After the conclusion of the LEP 
and Tevatron Run I experimental programs that reported no significant evidence for a Higgs 
boson, the experimental search for this particle has now become one of the major tasks of the 
Tevatron Run II and of the future LHC. Within the MSSM, the tree-level masses of the neutral 
Higgs bosons can be parameterized in terms of three input parameters: the mass of the CP-odd 
Higgs m A , the Z boson mass m z and the ratio of the two Higgs vacuum expectation values, 
tan/3 = v%jvx- At tree level, at least one of the MSSM Higgs bosons is bound to be lighter 
than the Z boson, thus the failure of detecting it at LEP indicates that the MSSM could be a 
realistic theory only after the radiative corrections to the Higgs boson masses have been taken 
into account. 

The radiative corrections arise from loop diagrams involving Standard Model particles and 
their superpartners. Although the earliest computations [3] of radiative corrections to the MSSM 
Higgs masses date back to the eighties, it was first realized in Ref. [I] that the inclusion of the one- 
loop top/stop 0{ott) corrections, where at = /if/(47r), h t being the superpotential top coupling, 
may push the light Higgs mass well above the tree-level bound. In the subsequent years, an 
impressive theoretical effort has been devoted to the precise determination of the MSSM Higgs 
masses: full one-loop computations have been provided leading logarithmic effects at two 

loops have been included via appropriate renormalization group equations El , and genuine 
two-loop corrections of 0(a t a s ) 0110111111121113, 0{a 2 t ) HHEEHHl, and 0(a b a s ) |E] have 
been evaluated in the limit of zero external momentum. The tadpole corrections needed to 
minimize the effective potential, V e g, have also been calculated ^01 at the same perturbative 
orders. Furthermore, the full two-loop corrections to the MSSM effective potential have been 
calculated in Ref. , together with a first study of the effect of the two-loop corrections to the 
Higgs masses controlled by the electroweak gauge couplings [T8] . 

The corrections controlled by the top Yukawa coupling are in general the most relevant ones. 
However, in regions of the MSSM parameter space where tan ^ > 1 the superpotential bottom 
coupling hb may be large (we recall that, at the classical level, hb/ht = (nib / nit) tan (5 ) and the 
one-loop bottom/sbottom corrections of 0(a\,), where a\, = h%/(4ir), can be numerically rele- 
vant and compete with those of Q(at). At the two-loop level, the evaluation of the corrections 
controlled by the bottom Yukawa coupling requires the inclusion of one-loop, tan /3-enhanced 
threshold corrections to the bottom mass JH]- If the physical bottom mass is used as input 
parameter in the one-loop part of the computation, potentially large tan /3-enhanced corrections 
appear at two loops. To address this problem, a set of renormalization prescriptions for the pa- 
rameters in the bottom/sbottom sector that avoid the occurrence of unphysically large threshold 
effects at two loops was proposed in Ref. ^5] for the 0(ai,a s ) part of the corrections. 

The purpose of this article is to complete the calculation of the two-loop, top/bottom Yukawa 
corrections to the Higgs boson masses in the effective potential approach. Such corrections were 
previously computed [HI I12| 11 4 j in the limit hb — > 0, which is accurate enough only for moderate 
values of tan/3. In that limit, the two-loop Yukawa corrections to the MSSM Higgs masses 
are of 0(a\m\), which we denote as 0(a\) for brevity. On the other hand, when the bottom 
Yukawa coupling is left arbitrary, the resulting two-loop corrections are proportional to various 
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combinations of couplings and masses: e.g., we find terms of 0(a 2 m 2 ), which might as well 
be interpreted as tan /3-enhanced terms of 0(atai,m 2 ). To simplify our notation, we will refer 
to all such "mixed" terms as to 0(at0b) corrections, and to the terms that depend only on the 
bottom Yukawa coupling as to 0(a 2 ) corrections. Our computation will thus provide us with the 
0(a 2 +atab+ot1) corrections to the MSSM Higgs masses, extending the 0{a 2 ) results obtained in 
Ref. [TI|. As a byproduct, we also calculate the 0(a 2 + atctb +ot 2 ) corrections to the minimization 
conditions of the effective potential. We express our results in the DR renormalization scheme, 
as well as in an "on-shell" scheme which extends the prescription described in Ref. ^3] to the 
case of the Yukawa corrections. The resulting analytical formulae are rather lengthy, thus we 
make them available, upon request , in the form of a Fortran code. 

The structure of this paper is the following. In section 2 we recall some general issues of 
the effective potential approach to the calculation of the Higgs masses. Section 3 describes our 
two-loop computation of the DR tadpoles and CP-odd, CP-even Higgs mass matrices, while 
section 4 addresses our on-shell renormalization prescription. Numerical results are given in 
section 5 and in section 6 we present a short discussion of the corrections controlled by the tau 
Yukawa coupling as well as our conclusions. 



2 Higgs masses in the effective potential approach 



We begin our discussion by recalling some general results concerning the computation of the 
MSSM Higgs masses in the effective potential approach. The effective potential, which we write 
from the start in terms of DR-renormalized fields and parameters, can be decomposed as V^fj- = 
Vq + AV, where Vq is the tree-level scalar potential and Ay contains the radiative corrections. 
Keeping only the dependence on the neutral Higgs fields H® and H®, the tree-level MSSM 
potential reads 



\H°\ 2 



\H° 2 \ 2 



(1) 



where: \i is the Higgs mass term in the superpotential (we assume it to be real, neglecting all 
possible CP-violating phases); irijj , TUjj and m| are soft supersymmetry-breaking masses; g 
and g' are the SU(2)l and £7(l)y gauge couplings, respectively. The neutral Higgs fields can 
be decomposed into their vacuum expectation values (VEVs) plus their CP-even and CP-odd 
fluctuations as H? = { Vi + Si + iPi)/V2 (i = 1, 2). The VEVs Vi are determined by solving the 
minimization conditions of the effective potential, i.e. 



dV cS 



dSi 



0, 



the second equality being automatically satisfied since we assume that CP is conserved. However, 



it is also possible to take v\ and v 2 as input parameters, or equivalently v z 



+ Vy and 



tan/3 = V2/V1, where v 2 is related to the squared running mass of the Z boson through m| = 
(g 2 +9 ,2 )v 2 /'i. In this case, the minimization conditions of V e s can be translated into conditions 
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on ijl 2 and mk: 



m 



17ln 



^ m 2 Hl + Si - (m 2 H2 + S 2 ) tan 2 (3 
2 + tan 2 /?-l 



— sin 2(3 +- tan 2/3 (mf^ 



m 



where the "tadpoles" Si and S 2 are defined as 

_ 1 <9AF 

Zj 7 ; = 



Vi dSi 



(3) 
(4) 

(5) 



In the effective potential approach, the mass matrices for the neutral CP-odd and CP-even Higgs 
bosons can be approximated by 



(Ml) 



eff 

ij 



d 2 V eS 



dPjdPj 



eff 

y 



d 2 V eS 



dSidSi 



(6) 



Exploiting the minimization conditions of the effective potential, Eq. P|). the CP-odd mass 
matrix can be written as 



M 



it 



-m 



2 V X V 2 



dPidPj 



(7) 



(A^p) e has a single non-vanishing eigenvalue that, in the approximation of zero external mo- 
mentum, can be identified with the squared physical mass of the A boson. We denote it as 
m 2 = m 2 + Am 2 , where m 2 = —2m 2 / sin 2(3 is the squared running mass of the A boson. The 
CP-even mass matrix can in turn be decomposed as 



(m 2 s ) CS = (M 2 S ) 



0,eff 



+ 



(am 2 ) 



eff 



(8) 



where the first term in the sum is the tree-level mass matrix expressed in terms of m z and m A : 

(M 2 )°' eS =( m z4 + m l s } -(m 2 z + m 2 A )spcp \ 
V s ) \ - (m 2 + m 2 ) s p cp m 2 s 2 p + m 2 A (% J ' 

{cp = cos (3 , sp = sin/3 and so on), while the second term contains the radiative corrections: 



(9) 



(AM 2 ) 



eff 



d 2 AV 



dSidSj 



d 2 AV 



dPidP, 



(10) 



It is clear from Eqs. (|7 |) -1|1U )) that, in order to make contact with the physical A mass, the 
effective potential should be computed as a function of both CP-even and CP-odd fields. 

Since V e fi generates one-particle-irreducible Green's functions at vanishing external momen- 
tum, it is clear that the effective potential approach neglects the momentum-dependent effects 
in the Higgs self-energies. The complete computation of the physical masses of the CP-even 
Higgs bosons, and m H , and of the CP-odd Higgs boson, m A , requires the full, momentum- 
dependent two-point functions. A detailed discussion of the correspondence between the effective 
potential approach and the full computation has been presented in Ref. Here we just notice 
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that the main conclusions presented in that paper regarding the 0(ata s ) and 0(af) corrections 
apply also to 0(a t af, + a^) corrections discussed here. Namely, the two-loop 0(0^0^ + 0^) correc- 
tions to the lightest Higgs eigenvalue are fully accounted for by the two-loop effective potential 
evaluation of supplemented by known momentum-dependent one-loop contributions, and the 
same is true for m H when m A is not too large. Instead, if m A > m t a full two-loop (D(atab + a\ ) 
computation of ra H requires additional momentum-dependent two-loop contributions, neglected 
by the effective potential calculation, that have not been computed so far. 

3 Computation of the two— loop Yukawa corrections 

We shall now describe our two-loop computation of the tadpoles Sj , the A-boson mass correction 
Am A and the matrix (A.M|) eff , including terms controlled by the top and/or the bottom Yukawa 
couplings. The computation is consistently performed in the gaugeless limit, i.e. by setting to zero 
all the gauge couplings, and by keeping ht and h\, as the only non-vanishing Yukawa couplings. 
In this limit, the tree-level (field-dependent) spectrum of the MSSM simplifies considerably: 
gauginos and Higgsinos do not mix; charged and neutral Higgsinos combine into Dirac spinors 
with degenerate mass eigenvalues |//| 2 ; the only massive SM fermions are the top and bottom 
quarks, while all other fermions and gauge bosons have vanishing masses; the only sfermions with 
non-vanishing couplings are the stop and sbottom squarks; the lightest CP-even Higgs boson, 

h, is massless, and the same is true for the Goldstone bosons; all the remaining Higgs states, 
(H, A, H^), have degenerate mass eigenvalues m? A . The tree-level mixing angle in the CP-even 
sector is just a = j3 — tt/2. 

The renormalization of the effective potential is performed according to the lines of Ref. |16| . 

i. e. we express V e fr, from the beginning, in terms of DR-renormalized fields and parameters. 
In practice, this amounts to dropping all the divergent terms in AV and replacing the two- 
loop integrals I(mf, m^m^) and J(m\,mi,) (see e. g. Ref. |16| for the definitions) with their 
"subtracted" counterparts I and J, first introduced in Ref. [20]. Alternatively, we could follow 
the procedure of Refs. |13U14| : express AV in terms of bare parameters and then renormalize the 
derivatives of AV (i. e. the tadpoles and the corrections to the Higgs masses), checking explicitly 
the cancellation of the divergent terms. The general formulae for the tadpoles and the corrections 
to the Higgs masses would look slightly more complicated in the latter case. However, we have 
checked that the two renormalization procedures lead to the same final result, as they should. 

According to Eqs. ©, Q and lfTU|) . the tadpoles and the corrections to the Higgs mass 
matrices can be computed by taking the derivatives of AV with respect to the CP-even and 
CP-odd fields, evaluated at the minimum of V c r. Following the strategy of Refs. ^JJ E], we 
compute AV in terms of a set of field-dependent parameters (masses and angles), and use 
the chain rule to express the corrections in terms of derivatives of AV with respect to those 
parameters. In each sector, the field-dependent parameters can be chosen as 

m q , m\ , m| , q , ip q , (p q , (q = t,b) , (11) 

where: m q and m~. are the quark and squark masses; 6 q is the field-dependent squark mixing 
angle, defined in such a way that < 6 q < tt/2 (to be contrasted with the usual field— independent 
mixing angle 9 q , such that — ir/2 < 9 q < vr/2); (p q is the phase in the complex quark mass; cp q 
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is the phase in the off-diagonal element of the squark mass matrix. For the explicit Higgs field 
dependence of these parameters, see Refs. ^JJE]- I* 1 the expression of AV relevant to the 0(a 2 ) 
corrections (i. e., with h b set to zero), the top and stop phases always combine in the difference 
ift — <pt , so that a convenient choice for the field-dependent parameter is c Vt -0 t = cos(c/?£ — (pt). 
On the other hand, when both ht and hb are nonzero, as it is the case in the present computation 
of the 0(a 2 +ottotb + a%) corrections, the situation becomes more complicated: besides the terms 
involving (ft — yp% and ipb — (fb, we find other terms, coming from diagrams with a charged Higgs 
or Goldstone boson, that involve the combinations ipt + (p b , fb + ft > ft + fb and fit + fb- 

Exploiting the field-dependence of the various masses and angles, we get the following general 
formulae for the 0(a 2 + ata b + ct b ) corrections in the DR renormalization scheme: 

off 1 

\ CA1 _ . O O — h _ . O . — h -L . o . O O — ,b 



(AM 2 s ) e u = 2hlmlF h l +2hlAbmbS 2 e b F h 2 +-h 2 bAlsl eb Fl 

+ ^h 2 fj 2 s 2 , et F^ + 2h t h b rnbfiS2e t F^ + h t hb^AbS2e t S2e b F 5 , (12) 

cff X 

(AMs) i2 = t$n m t s 2 e t F\ + -h 2 A t n s\ dt F\ + h t h b m b A t s 2dt Fl 

+ h 2 h iim h s 2 e b F\ + -hlAbU s\ db F| + h t h b m t A b s 2 e b F\ 

+ ^ h h s 2 8 t S28 b (A t A b + fi 2 ) F 5 + 2h t hbmtm b F 6 , (13) 

(AM 2 S )^ = 2h 2 t m 2 F t l +2h 2 A t mtS28 t F 2 t + ±h 2 A 2 t s 2 et Ft 

2 h l V 1 s le b F 3 +2h t h b m t n s 29b F\ + h t h b n A t s 2 e t s 2 e b F 5 , (14) 



+ 2 



vjfEi = mtn cot f3s 2 e t F t + m b A b S2e b F b + 2 m 2 b G b , (15) 



v 2 2 



S 2 = m b tan P s 2Sb F b + m t A t s 2St F t + 2 m 2 G l , (16) 



Ami = _J_ h ]»\ F t + f» A b F > + 2h t h b F A ) . (17) 
sn \ m~ — m~ mf — mf / 

p ^ \ h t% b\ 62 / 

In the equations above, At and A b are the soft supersymmetry-breaking trilinear couplings of 
the Higgs fields to the stop and sbottom squarks, and s 2 e q = sin 26g (q = t, b) refer to the usual 
field-independent squark mixing angles. The functions Ff (i = 1,2,3,4) , F 5 , F e , F q , G q and 
Fa are combinations of the derivatives of AV with respect to the field-dependent parameters, 
computed at the minimum of the effective potential; their definitions are given in the appendix. 
It can be noticed that, as it is predictable from the form of the MSSM Lagrangian, the above 
results are fully symmetric with respect to the simultaneous replacements t <-> b and Hi <-> H2, 
the latter resulting into tan/3 <-> cot/3 , v% <-> V2 , (AM 2 S ) C ^ <-> (AA^Dgf and Si <-> £2- 

An explicit expression of the two-loop top and bottom Yukawa contribution to AV can be 
found in Ref. [TJ], while the complete two-loop effective potential for the MSSM was given in 
the second paper of Ref. However, those expressions were computed for vanishing CP-odd 
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fields, thus omitting the dependence on the phases <p q and (p q . Since these phases appear in AV 
in many different combinations, it is not possible to infer the general field-dependent expression 
of AV by means of simple substitutions in Eq. (D.6) of Ref. ^2]) as it was the case in the 
computation of the 0(ata s + af ) corrections 2 . 

We worked out the general expression of the two-loop top and bottom Yukawa contribution to 
AV in terms of all the field-dependent parameters of Eq. 1)11(1 . We then computed its derivatives 
in order to obtain explicit formulae for the various functions appearing in Eqs. (|12j) - (|17|) . The 
use of a recursive relation for the derivatives of J(m| , m 2 ., rn^), presented in Ref. [TJi|, helped us 
to keep the number of terms involved under control. However, the resulting analytical formulae 
are very long and we choose not to display them in print. Instead, we make them available in 
the form of a Fortran code. 

4 On— shell renormalization scheme and input parameters 

The results presented in the previous section are valid when the MSSM input parameters are 
expressed in the DR renormalization scheme. This way of presenting the results is convenient for 
analyzing models that predict, via the MSSM renormalization group equations, the low-energy 
DR values of the parameters in terms of a set of boundary conditions assigned at some scale 
M GUT much larger than the weak scale (see Refs. |2H \'2'2\ for a list of public codes that are 
commonly used in this kind of analyses, and Ref. [23] for a comparison among them). General 
low-energy analyses of the MSSM, however, do not refer to boundary conditions at high scales, 
and are usually performed in terms of parameters with a more direct physical interpretation, such 
as pole masses and appropriately defined mixing angles in the squark sector. Such an approach 
requires modifications of our two-loop results, induced by the variation of the parameters entering 
the one-loop corrections when moving from the DR scheme to a different scheme (for a generic 
parameter x, we define the shift from the DR value x as 5x = x — x). 

However, it is not always possible to find a sensible definition with a direct physical interpre- 
tation for all the relevant parameters. For example, while there is a well known physical (= pole) 
definition for the masses, the so-called "On-Shell" (OS) definition, and an OS definition for the 
squark mixing angles can be also conceived [23], it is not clear what meaning should be assigned 
to an OS definition of parameters like {At , Af,, fi, tan (3). For instance, they could be related to 
specific physical amplitudes. However, given our present ignorance of any supersymmetric effect, 
such a choice does not seem particularly useful. In these cases it seems sometimes simpler to 
stick to a DR definition. 

It is rather easy to devise an OS renormalization scheme for the parameters in the top/stop 
sector, based on the OS prescription for the top and stop masses and the stop mixing angle 
and treating A t as a derived quantity, while retaining a DR definition for [i and tan/? (see 
e. g. Refs. I14j). Instead, some additional care is required in the choice of an OS scheme for 
the parameters in the bottom/sbottom sector, due to the potentially large one-loop threshold 
corrections ^5], proportional to tan/3, that contribute to the pole bottom mass. For example, 
a definition of A^ in terms of the OS bottom and sbottom masses and sbottom mixing angle, 
similar to the definition of At, would produce a shift SAf, proportional to tan 2 [3 [23]. When tan f3 
2 Also, we do not agree with Ref. j!2j on the sign of the penultimate line of Eq. (D.6). 
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is large, this would induce very large corrections to the Higgs masses at two loops, questioning 
the validity of the perturbative expansion. 

To overcome this problem, we adopt a set of renormalization prescriptions for the param- 
eters in the the bottom/sbottom sector, first introduced in Ref. ^S] for the case of the strong 
corrections, that avoid the occurrence of unphysically large threshold effects and at the same 
time enforce other desirable properties such as the decoupling of heavy particles, the infrared 
finiteness and gauge-independence. Generalizing these prescriptions to the case of the Yukawa 
corrections, and combining them with the usual prescriptions for the top/stop parameters 114] . 
we obtain a convenient OS renormalization scheme for the 0(atab + a b ) part of the corrections to 
the Higgs masses. Since the corrections controlled by the bottom Yukawa coupling can be size- 
able only for large values of tan (3, we work directly in the physically relevant limit of tan (3 — > oo, 
i. e. v\ — ► , V2 — ► v. 

For the OS squark masses and mixing angles, top quark mass and electroweak parameter 
v = (y^G^) -1 / 2 we adopt the definitions 

c 2 jrq / 2 \ ™ i nf 2 (mfj + nf 2 (m 2 - 2 ) y U^ w (0) 

(18) 

where q = (t,b), while IT?-(p 2 ) , £t(p) and Uy^ w (p 2 ) denote the real and finite parts of the 
self-energies of squarks, top quark and W boson, respectively. Following Ref. [Tl|, we further 
treat \i as a DR parameter computed at a reference scale Qo = 175 GeV, and ht and A t as 
derived quantities that can be computed by means of the tree-level formulae for mt and S2o t , 
respectively. In principle, we still have to define m , h and A b . However, in the large tan/3 
limit, the bottom mass is just zero, and the sbottom mixing angle becomes 

S28 b = — a - , 19 

m~ — m~ 
bi o 2 

which is independent of m and Af, . We can thus treat h as a quantity derived from the sbottom 
mixing, and use Eqs. (|T8|) and (fl9|) to obtain a prescription for 5hb'. 

Jft6 = M HzH + ^4_*l. (20) 

m bl ~ m b 2 S2e " 

In Ref. ^S], an OS definition for the quantity A\, = h A(„ or equivalently for 5Ab = {8A\, — 
5ht,Ai,)/hb , was proposed in terms of the (bib^A) proper vertex iAi2^(pf, p|, p\) for the case of 
the strong corrections. A generalization of that definition that can also encompass the Yukawa 
corrections is given by 

1 r ' ' -.2 r\\ I A / 2 



5A b = --j= \ K 12A (m~ bi , m- bi , 0) + A 12A (m^ , , 0) 



+ 



2 



H^m^-H^m^) H^m^-H^m^) UAAim^-UAAjm^) 

2 2 2 1 22 

m~ — m~ m~ — m~ m= — m~ 

Ol 62 0l 2 0l o 2 



(21) 
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Having fully specified our OS renormalization prescriptions in the limit tan/3 — ► oo, physically 
relevant for the 0(atctb + a l) corrections, we can proceed to obtain formulae for the CP-even 
Higgs masses in our OS scheme and merge them with the known 0(a 2 ) OS results |14j that 
contain an explicit dependence on tan/3. This can be done in three steps: first, we take the limit 
of tan/3 — ► oo , m — ► in the general DR results for the 0(a 2 +atab+a 2 ) part of (AAi 2 s ) e ; then 
we add the contributions due to the shifts of the parameters entering the one-loop corrections 
(this requires the computation of the 0(att + a ) part of the counterterms in the large tan/3 
limit); finally, we subtract from this results the pure 0(a 2 ) part which, being relevant for all 
values of tan f3, must instead be computed separately with the formulae of Ref. |14j . Notice that 
we do not encounter any terms that blow up when we take the limit of large tan (3 in the DR 
results: unphysically large contributions could only be introduced by hand, as the result of a 
poor choice of the renormalization conditions for the parameters in the bottom/sbottom sector. 

We discuss now the parameters that we will actually use as inputs of our calculation. In 
particular, although we have used Eqs. (|19j) - (f2~0~|) to define an OS bottom Yukawa coupling hf, 
through the sbottom mixing, we still need to exploit the experimental information on the bottom 
mass in order to obtain the DR running coupling hb. The OS coupling will then be computed 
through the relation hb = hb — <5/i . Following Ref. |15| . we define the running coupling hf, at the 
reference scale Qq = 175 GeV to be 

? _ > , n yDR _ g&vj / 99 ^ 
vi |l + e b | 

where: m = m ((5o)sM = 2.74 ± 0.05 GeV is the Standard Model bottom mass, evolved up 
to the scale Qq to take into account the resummation of the universal large QCD logarithms; 
e contains the tan /3-enhanced threshold corrections from both the gluino-sbottom and the 
higgsino-stop loops (denoted as e and e' h , respectively, in Eqs. (25) and (26) of Ref. [E]); 5b 
contains the residual threshold corrections that are not enhanced by tan/3. Notice that, as 
shown in Ref. }26| , keeping €b in the denominator of Eq. (|22[) allows to resum the tan /3-enhanced 
threshold corrections to all orders in the perturbative expansion. On the other hand, there is 
no preferred way of including the threshold corrections parametrized by 5b, whose effect on the 
value of hb is anyway very small. Neglecting all the terms controlled by the electroweak gauge 
couplings, e reads 



2 a s rrig \i tan (3 

q o 2~~ 

6 7r mf — mf 

01 02 



2 2 2 2 

mf mf mf mf 

bl ln-l-- B b2 s In 62 



at A t \i tan (3 
4 7r m~ — mk 

tl t2 



2 2 2 2 2 2 

mf — mf mf mf — mf mf 
bi y y 02 y y 

2 2 2 2 ~i 

mf mf mf mf 

*i In _ *L l n *a 



m| — fi 2 fJ? m~^ — fi 2 fi 2 



(23) 



It appears from Eq. ((2*2*]) that hb blows up when e approaches —1, in which case the correct 
value of the bottom mass cannot be reproduced with hb in the perturbative regime, and the 
corresponding set of MSSM parameters must be discarded. It can also be noticed from Eq. (|23|) 
that, since we take mg > , for At > (At < 0) the 0(a s ) and 0(at) contributions enter e with 
the same (the opposite) sign. Moreover, if we take only the 0(a s ) part of e into account, h\ 
can be enhanced by the threshold correction only for large values of m g and large and positive 3 
3 Our convention for the sign of fi is such that, e. g., the sbottom mixing parameter reads Xb — Ab + /i tan/3. 
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values of (a, whereas, when we include the 0(at) part, hb can be enhanced also for small values 
of rrig and large and negative values of both \i and A t . 

For the top/stop sector, we take as input the current central value of the top pole mass, 
rat = 174.3 GeV [27|, an d the parameters (toq i , m\j , A t ) that can be derived by rotating the 
diagonal matrix of the OS stop masses by the angle defined in Eq. ()18j) . Concerning the 
sbottom sector, additional care is required, because of our non-trivial definition of hf, and of the 
fact that, at one loop, the parameter m,Q ^ entering the sbottom mass matrix differs from the 
corresponding stop parameter txiq i by a finite shift We start by computing the renormalized 
coupling as given by Eq. ()20|) and ()22|) . Then we compute m,Q ^ following the prescription 
of Ref. • Finally, we use the parameters hb and ruq ^ to compute the actual values of the 
OS sbottom masses and mixing angle. Concerning the A-boson mass, which enters the tree- 
level mass matrix for the CP-even Higgses, we take as input the physical mass m A , dropping 
the distinction between m A and the effective potential mass fn A (this amounts to neglecting 
the effect of the uncomputed momentum-dependent two-loop corrections). The renormalization 
of the Z-boson mass, whose numerical value we take equal to the physical mass m z = 91.187 
GeV, does not affect the 0{atoib + a%) corrections. The remaining numerical inputs are the 
OS electroweak parameter v = 246.218 GeV and the strong coupling constant, which we fix as 
a s (Q ) = 0.108. 

5 Numerical results 

We are now ready to discuss the numerical effect of our two-loop corrections. In the previous 
sections we have discussed how to express our results in either the DR renormalization scheme or 
an OS scheme suitably chosen to separate the genuine two-loop corrections from the threshold 
corrections in the relation between mb and hb- Our DR results for the two-loop corrections to the 
Higgs masses and to the electroweak symmetry breaking conditions can be easily implemented 
in the existing codes [22] that compute the MSSM mass spectrum from a set of unified 
parameters at the scale M GUT . A study of the implications of our results in the framework 
of gravity (mSUGRA), gauge (GMSB) or anomaly (AMSB) mediated supersymmetry breaking 
models goes beyond the scope of this paper, and will appear elsewhere |28| . In the following 
discussion, we will adopt a low-energy point of view and assume that the various input parameters 
are related, when possible, to physical quantities. To this aim, we will make use of the OS 
renormalization scheme presented in section 4. We recall that, although our OS prescription is 
defined in the limit tan/3 — > oo, the corrections have an indirect dependence on tan/3 coming 
from the input value for hb , see Eq. (|22|) . 

In Figs.^, (upper panel) and^ (lower panel) we show the light Higgs mass as a function 
of tan/3, for m A = 120 GeV and m A = 250 GeV, respectively. The other input parameters are 
chosen as A t = 1 TeV, Aj, = 2 TeV, /i = tuq i = mu = rri£, = m g = 1 TeV. In each plot, the long- 
dashed curve corresponds to the value of obtained at 0(at + ata s + a^), i. e. by including only 
the one- and two-loop corrections controlled by the top Yukawa coupling; the dot-dashed curve 
includes in addition the one-loop O(o.b) corrections, controlled by the bottom Yukawa coupling 4 ; 

4 In the calculation of the one-loop O(at) and 0(af,) corrections we include the effects proportional to m% and 
the momentum corrections as in 
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Figure 1: The mass rrih as a function of tan/3, for m A = 120 GeV (upper panel) or 250 GeV 
(lower panel). The other input parameters are A t = 1 TeV, = 2 TeV, \i = mq i = mu = 
mo = rrig = 1 TeV. The meaning of the different curves is explained in the text. 
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Figure 2: The mass mh as a function of tan/3, for m A = 120 GeV, \i = At = —2 TeV, Ab = —3 
TeV, rriqi = my = tod = 1 TeV and rrig = 200 GeV. The meaning of the different curves is 
explained in the text. 

the short-dashed curve includes the two-loop 0{aba s ) corrections computed in Ref. ^5]; finally, 
the solid curve corresponds to the full two-loop Yukawa computation of m^, i. e. it includes also 
the O(ato.b + a%) corrections discussed in the previous sections. We can see from Figs.^i and^> 
that the corrections controlled by the top Yukawa coupling depend very weakly on tan/3 when 
the latter is large. On the other hand, the 0{czb) corrections lower considerably nih when tan/3 
increases. Concerning the two-loop corrections controlled by the bottom Yukawa coupling, the 
comparison between the dot-dashed and short-dashed curves shows that the 0(otba s ) corrections 
amount to a small fraction of the 0(otb) ones, but they can still lower by several GeV when 
tan (3 is large. The comparison between the short-dashed and solid curves shows that the effect 
of the 0{atoib + a?) corrections can also amount to several GeV when tan/3 is large. From 
Fig. we see that, when m A is small and the correction to rrih is mainly driven by (AA^g)n, 
the 0(atab + af) corrections enter with the same sign as the 0(a>ba s ) corrections, reducing 
further the value of m/j. On the other hand, Fig. ^> shows that for larger values of m A , when 
the correction to rrih is sensitive to (A«Mo)22j the new corrections account for an increase in nih 
of a few GeV at moderately large values of tan/3 (i. e., tan/? « 30-40). This is basically due to 
a positive contribution to (AA4|>)22 coming from the (D(atotb) part of corrections. When tan/3 
takes on larger values, however, the overall effect of the O(ato.b + a^) corrections to nih turns 
again to negative. 

It is interesting to realize that the 0{atab + a b 1 ) corrections can be sizeable also for parameter 
choices that make the 0(ctbCt s ) corrections irrelevant. In Fig. [21 we show as a function of 
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Figure 3: The masses and m H as a function of m A , for tan (5 = 40, At = 1 TeV, Ab = 2 TeV, 
[i = niq i = mu = rri£) = m g = 1 TeV. The meaning of the different curves is explained in the 
text. 

tan/3, for m A = 120 GeV, ^ = A t = —2 TeV, A\> = —3 TeV, trig j = mp = mo = 1 TeV and 
nig = 200 GeV. The meaning of the various curves is the same as in Fig. ^ Due to the small 
value of the gluino mass with respect to the sbottom masses, the 0(at,a s ) corrections to mh 
are negligible (in fact, the dot-dashed and short-dashed curves overlap). On the other hand, 
comparing the short-dashed and solid curves we see that the O(o.totb + a%) corrections can still 
amount to a few GeV when tan (5 is large enough. 

Finally, Fig. |3] shows both CP-even Higgs masses, and m H , as functions of the CP-odd 
Higgs mass, in the region of relatively small m A (100 GeV < m A < 180 GeV), for tan (5 = 40. The 
other input parameters are chosen as A t = 1 TeV, Af, = 2 TeV \x = tuq £ = mu = mo = mg = 1 
TeV. The meaning of the various curves is the same as in Fig.^ Comparing the short-dashed and 
solid curves we see that, for this choice of parameters, the effect of the 0(atctb + ay*) corrections 
is particularly evident in the region where m A is small, and can account for variations of several 
GeV (around 5 in this example) in both mh and m H . 

In general, it appears from Figs.Etllthat the two-loop 0(aba s ) and 0{atab + a b 1 ) corrections 
are usually a small fraction of the one-loop 0(ctb) ones. We stress that this is a desirable 
consequence of our renormalization prescription, which allows to set apart the tan /3-enhanced 
threshold corrections, resummed to all orders in the renormalized coupling fib- If we were to 
adopt for the bottom/sbottom sector the same renormalization prescription that we use for the 
top/stop sector, the dependence on tan/3 of the one-loop corrections would be smoother, but 
very large corrections would appear at two loops, questioning the validity of the perturbative 
expansion. 
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To conclude this section, we notice that our knowledge of the general formulae for the correc- 
tions to the CP-even Higgs mass matrix in the DR scheme allows us to estimate the uncertainty 
connected with the fact that we take the limit of tan (3 — * oo , mj, —* in the corresponding 

05 results. In the numerical examples considered above we find that, in the regions where the 
corrections are sizeable, the DR results for the 0(aba s ) part of (AM 2 S ) vary by less than 20% 
when the limit tan (3 — ► oo , mj, — > is taken. The DR results for the 0(atctb + a%) part of the 
corrections vary instead by less than 10%. We can assume that similar variations occurr also in 
the corresponding OS results, which leads to shifts in rrih typically smaller than 1 GeV. 

6 Conclusions and discussion 

In this paper we computed the 0{a. t ab + a%) corrections to the MSSM neutral Higgs boson 
masses and to the minimization conditions of the MSSM effective potential. Such corrections 
are relevant when the ratio of the Higgs VEVs, tan/3, is large. Combined with the previously 
computed 0{ata s ) [T^], O(a^) ^1] and 0(aba s ) 15 corrections to the neutral Higgs masses, 
and with the corresponding corrections to the minimization conditions of the effective potential 
[Tfij . these results provide us with a complete computation of the leading two-loop corrections 
controlled by the top and bottom Yukawa couplings. 

Using the formalism of the effective potential, we obtained complete analytic expressions for 
the momentum-independent part of the 0(af + atatb + af) corrections, valid for arbitrary values 
of the MSSM input parameters when the latter are expressed in the DR renormalization scheme. 
We also discussed an OS renormalization prescription for the parameters of the bottom/sbottom 
sector in the physically relevant limit of large tan (3. Such prescription, first introduced in Ref. |15j 
for the case of the 0(aba s ) corrections, allows to separate the large threshold corrections ap- 
pearing in the relation between hb and the pole bottom mass from the genuine two-loop effects. 
Finally, we discussed the numerical impact of our results in a few representative examples, show- 
ing that, for large values of tan/3, the 0{atdb + af) corrections can induce variations in the 
Higgs masses of the order of a few GeV. Since our analytic expressions, both in the DR and the 
OS schemes, are too long to be useful if explicitly written on paper, we choose to make them 
available in the form of a Fortran code. 

Although the terms controlled by the top and bottom Yukawa couplings undoubtedly account 
for the bulk of the two-loop corrections, several pieces are still missing for a complete two- 
loop computation of the MSSM Higgs masses. When tan (3 is large, the corrections controlled 
by the tau Yukawa coupling h T might in principle be non-negligible. In the approximation 
of neglecting the electroweak gauge couplings, the only two-loop corrections involving the tau 
Yukawa coupling are those of (D(aba T ) and those of 0(a?), where a T = h%./{4nv). While the 
mixed 0(ctba T ) corrections would require a dedicated computation 5 , explicit formulae for the 
0{oc^) corrections can be obtained from the formulae of Refs. (141 I16j for the purely O(a^) 
corrections, with the replacements t — ► r , bi — > v T , N c — > 1 and H\ <-> Hi [the latter resulting 
into tan/3 «-> cot (3 , v\ <-> i>2 , (AM S )°^ (AM. S ) C ^ and Si «-> £2]. If the input parameters 
are given in the OS scheme, a suitable definition of 5A T is required in order to avoid introducing 
tan /^-enhanced terms in the two-loop part of the result. Anyway, we find that the 0(a^) 
5 We thank A. Brignole for drawing our attention on this point. 
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corrections to the Higgs masses are in general very small compared with those controlled by the 
bottom Yukawa coupling. Besides the hierarchy between and m T , the suppression of the tau 
corrections is motivated by the absence of color enhancements, and by the fact that the only 
tan /^-enhanced threshold corrections to the relation between h T and m T are those controlled by 
the electroweak gauge couplings. 

A full two-loop determination of the MSSM Higgs masses will require going beyond the 
gaugeless limit and the effective potential approximation, i. e. including both the corrections 
controlled by the electroweak gauge couplings and the effect of the momentum-dependent part 
of the Higgs self-energies. It can also be expected that, among the three-loop corrections, at 
least those involving the top Yukawa coupling affect the Higgs masses in a non-negligible way. In 
Ref. ^S] the two-loop, zero-momentum electroweak corrections have been computed numerically 
in a representative scenario, and found to yield a shift in the lightest Higgs boson mass of 
about 1 GeV with respect to the result obtained in the gaugeless approximation. In Ref. j2H], the 
theoretical uncertainty in the prediction for rrih arising from the combined effect of the missing 
two-loop corrections and the leading three-loop corrections has been estimated to be around 3 



If the MSSM is a viable theory for physics at the weak scale, a light Higgs boson will be 
discovered either at the Tevatron or at the LHC. Subsequently, its properties will be determined 
with high precision at a future linear collider: for example, the predicted experimental accuracy 
in the determination of m,h at TESLA is about 50 MeV [HO] • It is thus clear that further effort will 
be required in the coming years, in order to improve the accuracy of the theoretical predictions 
up to the level required to compare with the experimental results expected at the next generation 
of colliders. 
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Appendix 

We present here the expressions for the functions Ff (i = 1,2,3,4) , F$ , Fq , F l , G l and Fa, 
appearing in Eqs. (|12)) - (|17|) . in terms of derivatives of the DR-renormalized AV^, computed at 
the minimum of V e s: 
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In the above formulae, z q = sign(A g ), 6 where X q (q = t, b) is the squark mixing parameter. The 
functions Ff , F b and G b can be obtained from their top counterparts through the replacement 
t <r-> b. 
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